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The de Broglie-Bohm approach permits to assign weU defined trajectories to particles that obey 
the Schroedinger equation. We extend this approach to electron pairs in a superconductor. In the 
stationary regime this extension is completely natural; in the general case additional postulates are 
required. This approach gives enlightening views for the absence of Hall effect in the stationary 
regime and for the formation of permanent currents. 
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I. INTRODUCTION 



If the evolution of a particle obeys the Schroedinger equation, then the de Broglie-Bohm quantum theory (dBB)i*2ii 
provides a deterministic description of its motion by assuming that, besides the classical forces, an additional "quantum 
force" V(fi.^V^|'0|/2m|'0|) acts on the particle. The dBB theory remains applicable in the presence of magnetic fieldsi^ 
Variations and extensions of dBB have recently been considered f^i?. 

In this article we want to extend dBB to the motion of Cooper pairs in a superconductor, which obey the Ginzburg- 
Landau equation2i^ 

^(-iW-|A)% + a7^ + /?|V'|'^ = -7?i^ • (1) 

Here we follow the notation of Ref. |^ that makes Eq. Q look similar to the Schroedinger equation: ip is the "order 
parameter" , t is the time, q and m are the charge and the mass of a pair, h and c are Plank's constant and the speed 
of light, a, (3 and 7 are material constants, and A is the vector potential of the electromagnetic field. We have chosen 
a gauge such that the scalar potential vanishes. 

In this letter we will interpret the order parameter 'ip — {iple'^'^ as the vawefunction for Cooper pairs. As usual, 
will have the meaning of density of particles, h^ip will be the canonic momentum, and v = {h^ip — 2A)/m, the 
velocity of the pairs. However, the equation of motion will not be based on the usual Hamilton-Jacobi formalism, 
but just on Newton's second law. The price we will have to pay for this approach is some basic knowledge of fluid 
mechanics (e.g. Ref. |9|). 

II. STATIONARY SITUATIONS 

In this case the right hand side of Eq. vanishes. The remaining equation differs qualitatively from the 
Schroedinger equation, since it is nonlinear. We shall see that, nevertheless, the dBB theory can be smoothly 
extended. 

A. Equation of Motion 

For a stationary situation, conservation of particles gives V • (jV'I'^v) ~ 0, which leads to 

2{hVip ~ -A) • VIV'I + iV'KW^vs - -V • A) . (2) 

c c 

Expanding Eq. using Eq. jSJ, and defining 
we arrive at 

-Qstat = {hV^ - ^ A)2 + a . (4) 
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Taking the gradient at both sides of Eq. and noting that in the stationary regime a = (v • V)v is the convective 
time derivative of the velocity, and hence the acceleration of a pair, we obtain 

V 

- VQstat + g- X B = ma , (5) 

c 

with B = V X A. This means that, in addition to the magnetic force, the quantum force —VQstat acts on every 
pair. As a matter of principle, —VQstat is not only the quantum force, but it also includes the force exerted by the 
lattice. Note that the only physical input required for Eq. JSJ is Eq. (with vanishing r.h.s.) and the interpretations 
made in the Introduction; the material parameters, the applied field or the geometry of the superconductor have no 
influence. 



B. The Hall Effect 



We shall see that Eq. ||SJ) provides an intuitive way to explain the absence of Hall effect in a superconductor 
in a stationary state. The Hall effect was discovered in 1879, well before quantum effects were suspected to exist. 
Historically, it turned out that an apparently trivial classical explanation could be given to the effect (except for the 
existence of positive carriers). Only when the Hall voltage turned out not to be a linear function of the magnetic 
fieldjiS the effect was called "quantum" . 

In the case of a superconductor carrying a stationary current perpendicular to a magnetic field that penetrates the 
sample, there is clearly a Lorentz force that acts on the flowing electrons; therefore, it seems obvious to anticipate 
the presence of an electric force that balances it, implying a Hall voltage. 

However, early experiments found that the Hall voltage drops to zero when a sample becomes superconducting^Sii^ 
Moreover, if there were an electric field in the superconductor, it would produce normal currents and dissipation. At 
least in the case of permanent currents, this is not what is observed. 

The resolution of the paradox is provided by Eq. jSJ): the Lorenz force is balanced by the quantum force —VQstat 
in order to complete the total force required to keep the pairs on their trajectories, and no electric force is required. 

Later experiments did find a Hall voltage in superconductors (e.g. Ref. [ij)- However, in these cases the voltage is 
due to the motion of vortices and cannot be described as a stationary situation. 



C. An Example 

As an illusration, let us examine a situation with constant field and cylindric symmetry. Let the field be B = Bqz 
and the vector potential, A = ^Bord, where r, 6 and z are cylindric coordinates. Cylindric symmetry implies that (p 
is a linear function of 9. Together with single- valuedness, this enables us to write tp — —n9, where n is an integer. It 
follows that V — {l/m){\q\Bor/2c — nh/r)9 and therefore the sum of the Lorentz and the centrifugal force is 

_ I (n^h^ q^B^r"^ 

" Lorentz ^ centrifugal 



mr 



4c2 



This force points away from or towards the z-axis, depending on whether r is smaller or greater than rs — 
{2nch/\q\Boy/^. On the other hand, using Eq. we find 

_ 1 fnh \q\BQr^^ 

— ttfstat — -Z 



2m \ r 2c 

(The quantum potential is maximal at r = rs-) We can now immediately verify that —VQstat = — FLorcntz 

^centrifugal- 



III. GENERAL CASE 



Now Eq. ^ is not only nonlinear; it is also nonunitary. Accordingly, the number of Cooper pairs is not conserved 
and any model that is limited to a fixed number of par ticles cannot describe the physical situation. We can still 
envision a deterministic evolution in the spirit of R.efs. \lH andlTfil ( "beables" ). The objects of the model are the order 
parameter ip^ which acts as a field, and the density of particles. We still postulate that particles move with velocity 
(h^if — ^A)/to. The order parameter, the density and the velocity depend on position and on time. The evolution 
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of the system will consist of two processes: motion of pairs and creation (or destruction) of pairs. In order to have 
a consistent model, we require that if the density of pairs is initially IV'Pi it will remain given by ji/jp as the system 
evolves. 



A. Pair Trajectories 

We expand Eq. ^ and factor out e*"^. Equating the imaginary parts we obtain 



V-(|^pv)^27|^P^ (6) 



and, from the real parts 



We define now 



Taking the gradient at both sides of Eq. {Tj), using Eq. ©, noting that the acceleration is now a = (v • V)v + dv/dt 
and the electric field is E = —{l/c)dA/dt gives 



-V(Qstat + (9dyn)+g(E+- xBj =ma, (9) 

so that the motion of the pairs can still be interpreted as having Newtonian trajectories, provided that the quantum 
force — V(Qstat + Qdyn) IS added to the electromagnetic force. Note that now the quantum force is not a function of 
lipl only, but is still a functional of '0- 

In regions where \^\ is nearly constant and uniform, — VQdyn ~ this case — VQdyn behaves as a viscous 

force, with viscosity coefficient ^It/i^. 



B. Creation and Annihilation of Particles 



We build a model for formation and destruction of pairs as follows. We define the depairing potential as 

Qdcp = 2 (^7 + (g^tat +a+ y w^) + 27gdyn ■ (10) 

We postulate that at intervals of time r the system checks itself. (For typical non-stationary processes in superconduc- 
tors, characteristic times arc of the order of 10~^^s.) At the end of each interval we divide the sample into "depairing 
cells" such that the volume integral in each cell obeys | / Qdcp\'4'\'^dV\ = H/t. Since the integral of QdcplV'l^ over the 
entire sample is not necessarily integer, some marginal region of the sample will in general be left out of the cells. In 
order to have a well defined division we need some additional criterion; for instance, we may require that the average 
boundary area of the cells be as small as possible. This criterion would discourage elongated cells or merging of cells 
with positive and negative Qdcp; also, the regions left out of the cells will be those where IQdopHV'l^ is small. Finally, 
we postulate that if / Qdcp\ip\'^dV is negative, then a new Cooper pair forms at the center of the cell and, if it is 
positive, the pair closest to the center is destroyed. 

We have to prove that the evolution of the pair density predicted by this model is the same as that predicted by 
Eq. JQ). The pair density changes due to two processes: flow of pairs and net pair creation. The increase of density 
per unit time do to flow is —V • (jV'I'^v). The increase due to creation equals the net number of negative cells per 
unit volume, — QdopIV'!'^'''/^! multiplied by the frequency 1/r. (We have assumed that the cells and r are sufficiently 
small to be treated as a continuum, with QdoplV'P uniform in the analyzed region.) Substituting Eqs. {T)) and © 
into Eq. l{TU\i we obtain that the total rate of density increase is d\tp\'^/dt, as required. 

For fast processes we expect that typical depairing cells will have microscopic sizes, but if the particle density is 
small or a process is almost stationary, these cells can be large. We could build an alternative model by postulating 
the existence of depairing cells with fixed positions and sizes, which destroy or create a pair every time that the 
integral / / Q dcp\'4'\'^ dV dt in the cell changes by ±h. 
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C. Example 

We consider ring of radius R, sufficiently thin to be treated as one-dimensional, threaded by a magnetic flux <i>. For 
negative times the ring is in the normal state, but at t = it is instantaneously cooled and becomes superconducting. 
In Eq. cooling below the critical temperature is implemented by changing the value of a from positive to negative. 
It is well known that if $ is not an integer multiple of the quantum of flux $0 = ch/\q\, then a permanent current 
flows around the ring. A classical-minded question naturally arises: since there are no electromagnetic fields in the 
ring, how do the charges know to start moving when the temperature is lowered? 

In order to answer this question we first have to calculate tp- Let us restrict our attention to the period during which 
\tp\ is so small that the nonlinear term in Eq. can be neglected. For a thin ring Eq. takes the one-dimensional 
form 

.d $ V , , ^ dij 

I— ip + aip = —nj—- . 



We write ■0 = J2'^=-oo Tn{t)e^'^'^^ and linearity permits to treat each harmonic separately. We obtain 



Since the ring was initially in the normal state, the initial values T„(0) are very small. (They don't vanish due to 
thermal fluctuations.^^) 

For simplicity, let us consider the case in which An > for just one value of n (the closest to $/$o), which will 
be denoted by n. This means that Tn will increase with time, whereas r„ will remain negligible for n ^ fi. As a 
consequence, [V'l will be independent of position. 

After Tp is known, we can evaluate the potentials. We obtain Qstat = 0, Qdyn = fin^a and Qdep — —2hXn < 0. 
Since these potentials are independent of position, it follows that the pairs are not accelerated. Current arises due to 
their formation, with velocity {h/mR){^/^o — n)9. 



IV. CONCLUSION 



We have extended the de Broglie-Bohm theory, which was built for particles that obey the Schroedinger equation, 
to the case of electron pairs in a superconductor, which obey a more complex (nonlinear and nonunitary) equation. 
For the stationary regime this extension is completely natural; in the general case, in which the number of pairs is 
not conserved, additional postulates are required. 
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